Soliton Solutions of the KP Equation with V-Shape Initial Waves 
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We consider the initial value problems of the Kadomtsev-Petviashvili (KP) equation for symmetric 
V-shape initial waves consisting of two semi-infinite line solitons with the same amplitude. Numerical 
simulations show that the solutions of the initial value problem approach asymptotically to certain 
exact solutions of the KP equation found recently in [1]. We then use a chord diagram to explain 
the asymptotic result. We also demonstrate a real experiment of shallow water wave which may 
represent the solution discussed in this Letter. 
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We consider the Kadomtsev-Petviashvili equation in 
the form, 

(4ut -I- Quu^ + Uxxx)x + 3u.y.y = 0, (1) 

which describes quasi-two dimensional, weakly nonlin- 
ear and long waves such as shallow water waves and ion 
acoustic waves PJ |3] 0]. It is well-known that the KP 
equation admits iV-soliton solutions, each of which has 
asymptotically the same set of N line-solitons in both 
y ±oo. Interaction property of two line-solitons with 
the same amplitude has been studied based on a two- 
soliton solution which forms an "X" shape with a phase 
shift due to the interaction. We call this two-soliton so- 
lution 0-type, where "O" stands for original. However, 
in 1977, Miles [5| pointed out that 0-type soliton so- 
lution becomes singular if the angle of the interaction is 
smaller than certain critical value. Since the KP equation 
is supposed to give better approximations in this regime 
of small angles, it sounds strange that we do not have a 
reasonable solution of the equation. Miles also found that 
at the critical angle the two line-solitons of the 0-type 
solution resonantly interact, and a third wave (soliton) is 
created to make a "Y" -shape solution, which is indeed an 
exact solution of the KP equation. After the discovery 
of the resonant phenomena in the KP equation, several 
numerical and experimental studies were performed (see 
e.g. [SI III E]). However, no significant progress has 
been made in the study of the soliton solutions of the 
KP equation for almost quarter century. 

In the last five years, a large variety of new soliton so- 
lutions has been found and classified [1] [101 HH IHl US] ■ 
Those new soliton solutions enable us to describe the 
interaction properties of line-solitons even in the param- 
eter regime of small angles where 0-type solution be- 
comes singular. In this Letter, we report how some of 
these new solutions appear under certain physical set- 
tings considered in the studies on the generation of freak 
(or rogue) waves in shallow water [HI [151 [IS]- We also 
present an elementary experiment of shallow water wave 



demonstrating a real existence of those new solutions. 

Let us first recall that u(x,y,t) = 2\hi T{x,y,t)]xx is a 
solution of ([T]) , if the r-function r is given by the Wron- 
skian form, 

T{x,y,t)^Wr{h,--- Jn), (2) 

where the functions /^'s satisfy the /mear equations, fy — 
fxx, ft = -fxxx (see e.g. [TT]). The soliton solutions we 
consider are generated by 

Here the coefficient matrix A — (a^j) is a constant TV x 
M matrix, and kj are constants with the ordering ki < 
■■■ < kf^j. Thus each solution u{x,y,t) is completely 
determined by the A-matrix and the /c-parameters. 

One soliton solution is then obtained in the case with 
iV = 1 and M = 2: We have r = e^' -I- e^^ for < kj, 
and u — 2(lnT)a;x gives 

u^^-{k,-hfsech'^-i0,-e,). (4) 

The line 6i — 9j represents the ridge of the soliton, and 
we refer to this Zme-soliton as [i, j]-soliton. The ampli- 
tude Q:[i, j] and the inclination 7[j, j] := tan?/;[i, j] of this 
soliton are given by 

= \ikj - hf, 7[i, j] = fcj -f kj, (5) 

where — 7r/2 < j] < Tr/2 is the angle of the line soliton 
measured counterclockwise from the y-axis. This angle 
also represents the propagation direction of the line soli- 
ton (see FigjT|. 

In the case with iV = 2,M = 4, using the Binet- 
Cauchy theorem for the determinant, the r-function 
r = Wr(/i, /2) can be expressed in the form, 

r= J2 (6) 

l<i<j<4 
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FIG. 1: Contour plots and the chord diagrams of sohton so- 
lutions. Left: O-type solution. Right: (3142)-type solution. 
Each [i, ji] denotes the [i, j]-soliton. The length of [1, 4]-soliton 
changes in t. The upper (lower) chords represent the solitons 
in y > (y < 0). The thicker chords correspond to the soli- 
tons in the right side (a; > 0). 



where is the 2x2 minor consisting of i-th and j-th 

columns of the matrix A, and E{i,j) = {kj — ki)e^''^^^ . 
For the regular solutions, we require all of these minors 
to be non-negative (note E{i,j) > 0). 

As was shown in [TJ [T2J [T^ , each r-function ^ gen- 
erates a soliton solution which consists of at most two 
line-solitons for both y ±oo. In the cases where the 
T-function ^ has at least 4 terms (i.e. A is irreducible, 
see [1, 12 ), we have seven different types of soliton so- 
lutions. Two of them are usual two-soliton solutions, 
O-type and P-type (this type fits better in the physical 
setting for the KP equation), and they are steady trans- 
lational solutions with an "X" -shape. The other five are 
non-stationary solutions. 

In particular, we consider the following two types 
which are relevant to the solutions of the initial value 
problems considered in this Letter: One is the O-type 
solution which consists of two line-solitons of [1,2] and 
[3,4] for y — > ±oo. Other one is non-statinary, and it 
consists of [1,3] and [3,4] line-solitons for y +oo and 
[1,2] and [2,4] line-solitons for y — oo. Let us call this 
soliton (3142)-type, since those four line-solitons repre- 
sent a permutation vr = (3142)- Figure [I] illustrates 
the contour plots of O-type and (3142)-type solutions in 
the a;?;-plane, and the corresponding chord diagrams rep- 
resenting each soliton as a chord joining a pair of fc^'s. 
Note in particular that the length and the location of 
each chord give the amplitude and the angle of the cor- 
responding soliton (cf. ([5])). The A- matrices for those 
solutions are respectively given by 



An^ 



1 a 
16 



A 



(3142) 



1 a 
1 



(7) 



where a, 6, c > are constants determining the locations 
of the solitons (see [13]). Notice that the r-function for 



the (3142)-type contains 5 exponential terms, 
r=(/!c3-fci)e''i+'^-'-|-(fc4-fci)6e''i+''^ + (fc3-fc2)ae^"+^' 



-I- ( fc4 - fca ) a6e''^ -I- ( fc4 - fca ) ce' 



3 + 04, 



(8) 



and the r-function for O-type with c = in ([8| contains 
only 4 terms. 

Let us fix the amplitudes and the angles of the solitons 
in the positive x regions for both O- and (3142)-types, 
so that those solutions are symmetric with respect to the 
cc-axis (see Fig|T]): 



a[l,2] = a[3,4] (for O - type) 
a[l,3] = a[2,4] (for (3142) - type) 

-7[1,2]=7[3,4] (for O-type) 
-7[1,3]=7[2,4] (for (3142) - type) 



(9) 
(10) 



Then from ([5|, one can find the A:-parameters in terms 
of a and 7 with fci = —k^ and ^2 = — ^3 (due to the 
symmetry): In the case of O-type, we have 



ki = -7/2 



k2 = -7/2 




The ordering k2 < k^ then implies 7 > \/2a. 

On the other hand, for the (3142)-type, we have 



fci 



-7/2 



7/2 



(11) 



(12) 



The ordering k2 < k^ implies 7 < y/2a. 

Thus, if all the solitons in the positive x -region have 
the same amplitude a for both O- and (3142)-types, then 
an O-type solution arises when 7 > ^/2a, and a (3142)- 
type when 7 < \/2a. Then the limiting value at k2 = k^ 
{— 0) defines the critical angle, i.e. 



7c 



2a. 



(13) 



Note from (|8| that at the critical angle, i.e. k2 = fca, 
the r-function has only three exponential terms, and this 
gives a "Y" -shape resonant solution as Miles noted ^5^. 

One can also see from ^ that for (3142)-type solu- 
tion, the solitons in the negative x-region are smaller 
than those in the positive region, i.e. a[3,4] = a[l,2] = 
7^/2 < a, and the angles of those in the negative x- 
regions do not depend on 7 and 7 [3, 4] = —7 [1,2] = 
7c. Two sets of three solitons {[1, 3], [1,4], [3,4]} and 
{[2,4], [1,4], [1,2]} are both in the soliton resonant state, 
that is, they are resonant triplets. These properties of 
the (3142)-type solution are the same as those of Miles' 
asymptotic solution for the Mach refiection of a shallow 
water soliton f5^. 

We now discuss the numerical simulations of the KP 
equation with V-shape initial waves with a = 2 and one 
free parameter 7 (0 < 7 < 7r/2), 



u{x, y,0) — 2 sech^(a; — j\y\)- 



(14) 



This initial wave form has been considered in the study 
on an oblique refiection of a solitary wave in two layer 
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FIG. 2: Completion of the chord diagrams. The left diagrams 
show the initial V-shape waves with the same amplitude for 
O-type (top) and (3142)-type (bottom). The right diagrams 
represent the asymptotic solutions. 
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FIG. 3: The contour plots of the solutions. Upper: Nu- 
merical solution for V-shape initial wave with a = 2 and 
7 = 1.5. Lower: (3142)-type solution at the same times, and 
the parameters are (fci, . . . , k^) = (—7/4, —1/4, 1/4, 7/4) and 
{a,b,c) — (4,4/7,4/3) in the ^-matrix, so that all solitons 
meet at the origin at i = (see [13j). 

fluid due to a rigid wall along the x-axis [THlfT5] . The nu- 
merical simulations are based on a spectral method with 
window-technique where the boundary of the computa- 
tion domain is patched by the corresponding one-soliton 
solutions (the details will be published elsewhere). 

In terms of the chord diagrams, the V-shape initial 
waves are represented by the partial chords in Fig[T]with 
the thicker chords. Then the main goal in this Letter is 
to show that the asymptotic solutions of the simulations 
are given by O-type solution for 7 > 7c and (3142)-type 
for 7 < 7c. Namely, the partial chords are asymptot- 
ically getting to be the corresponding complete chords 
(see Figj2] and the discussion in the end of this Letter). 

Figure |3] shows the contour plots of the numerical so- 
lutions for 7 = 1.5 (< 7c) and the (3142)-type solution 
with the same 7. Note that the waves in the trailing part 
of the solution are fully developed at t = 20 within the 
sight of numerical domain. At i = 20 the numerical so- 
lution is in remarkable agreement with the (3142)-type 
solution: The waves in the trailing part in the negative x- 
region are almost identical with [3,4]- and [1, 2]-solitons 
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FIG. 4: Left: The time evolution of the stem amplitude Ura 
(solid line) for a = 2, 7 = 1.5. The dashed line denotes 
a[l,4] = 6.125. Right: The time evolution of the stem length 
for a = 2 and 7 = 1.5, 1.748, 1.898 < 7c = 2. The dashed 
lines denote ( 16 1. 



in y > and j/ < respectively. From the corner of 
the V-shape, a new wave is generated. This wave is 
called stem, and can be identified as the [1, 4]-soliton 
(see Fig|T]). The time evolution of the stem amplitude 
Um is shown in Fig|4] and it approaches to the asymp- 
totic value, a[l,4] = (7 + ^/2af /2 = 6.125. The time 
evolutions of the stem length Lg are also shown in Fig|4] 
Identifying the stem as [l,4]-soliton, the formula of Ls 
can be derived. The ridges of [1,3] and [1,4] solitons are 
given by 9i = ^3 and 9i = 6*4, which lead to 



1 



72/ 



(7c' + 37')i, X 



1 



4^ ,c . - , . , 4 

Then the stem length Ls is given by L 



ilc + lft. (15) 
2y with y at 



the intersection point of those lines (15 1, i.e. 

Ls = (7c - l)t- 



(16) 



We performed simulations for the cases with 7 = 
1.5, 1.748, 1.898 < 7c, and found that the solutions 
agree generally well with the corresponding (3142)-type 
solutions. However, as seen in Fig|4] the stem length 
Ls of the numerical solution is slightly longer than the 



expected value (16 1 when 7 approaches to 7c. 

For 7 > 7c = 2, we performed the simulations for the 
cases with 7 = 2.1,2.207,2.367,2.5, and confirmed that 
the solutions approach to O-type solutions with the same 
7. Namely the waves in the trailing part in cc < are 
identified as [1,2]- and [3, 4] -solitons with a phase shift 
5x which is given by [5] 



4 = 7c-'M7V(7'-7f)]- 



(17) 



Also using the exact solutions of O-type and (3142)-type, 
one can obtain the formula Ua in terms of 7 [5j, 



(7 + 2)2/2 for 7<7c 
8/(1 + 6-'^-) for 7>7c 



(18) 



Figure [5] illustrates the numerical results for the phase 
shift and the asymptotic maximum amplitude Ua with 



the formulae (17 1 and (18) 



At the critical angle 7 ~ 
of u is given by Um = 4q; 



7c -2, 
= 8 (cf. 



the maximum value 
( 1 8 1 ) . The numeri- 
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FIG. 5: Left: The phase shift &^ for 0-type solution. The 
dashed hne denotes p7| . The circle denotes numerical results 
for 7 = 2.1, 2.207, 2.367, 2.5. Right: The asymptotic value 
of the maximum of u{x,y,t). The dashed line denotes (18 1. 
The circle denotes numerical results. 




10 
t 



100 



FIG. 6: Left: The time evolution of the maximum Um of 
u{x,y,t) at the critical angle 7 = 7c = 2. Right: The time 
evolution of the phase shift 5x at the critical angle 7 = 2. 



initial V-shape wave as a (sub) chord diagram consisting 
of two chords which correspond to the semi-infinite soli- 
tons of the V-shape. Then we observed that the asymp- 
totic solution is given by a chord diagram which is the 
completion of the sub chord diagram. Here the comple- 
tion means that the resulting chord diagram should give 
the unique permutation and the corresponding solution 
has the smallest total length of the chords. For exam- 
ple, for the (sub) chord diagram of V-shape with 7 > 7^, 
there are two complete chord diagrams having this sub 
diagram (they are of 0-type and (2413)-type). The total 
length of solitons in (2413)-type is larger than the 0-type 
one. Thus the completion of the sub chord diagram is 
the diagram of the 0-type. This completion can be also 
applied for more general cases of V-shape initial waves 
with different amplitudes. The details will be reported 
elsewhere. Finally, in Fig|7] we present an elementary 
(desktop) experiment of a shallow water wave. The size 
of the tank is 37cm x 65cm with 2cm depth. The waves 
are generated by shifting the tank in a diagonal direction. 
The resulting wave is similar to (3142)-type. The exper- 
iment is easy, and of course, we can also get "Y" -shape 
and 0-type as well. 
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FIG. 7: A real experiment showing a (3142)-type solution. 

cal simulation shows that the amplitude approaches very 
slowly to the asymptotic value 8 (Figjoj). In the limit 
7 ^ 2+ of 0-type solution, the phase shifts become in- 
finity. On the other hand, in the limit 7 — > 2~ of the 
(3142)-type solution, the stem length Lg approaches to 
zero. We could not find an exact description of the so- 
lution at the critical angle, but the simulation indicates 
that the phase shift (or the stem length) seems to have 
a logarithmic increase (Figj6|. 

Let us now make a summary and some discussion: We 
performed several numerical simulations of the KP equa- 
tion with symmetric V-shape initial waves. Then we 
found that the solutions are asymptotically approaching 
to either (3142)-type solutions or 0-type solutions de- 
pending on the angles of the V-shape. Those solutions 
are expressed by the chord diagrams, equivalently the 
permutations (see FigjlJ. As shown in ^13J, the soliton 
solutions generated by the r-functions in the form ^ are 
uniquely expressed by the permutations. Then in order 
to explain the asymptotic solutions, we first express each 
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